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In order to model random density-dependence in population dy- 
namics, we construct the random analogue of the well-known lo- 
gistic process in the branching process' framework. This density- 
dependence corresponds to intraspecific competition pressure, which 
is ubiquitous in ecology, and translates mathematically into a quadratic 
death rate. The logistic branching process, or LB-process, can thus 
be seen as (the mass of) a fragmentation process (corresponding to 
the branching mechanism) combined with constant coagulation rate 
(the death rate is proportional to the number of possible coalesc- 
ing pairs). In the continuous state-space setting, the LB-process is a 
time-changed (in Lamperti's fashion) Ornstein-Uhlenbeck type pro- 
cess. We obtain similar results for both constructions: when natural 
deaths do not occur, the LB-process converges to a specified distri- 
bution; otherwise, it goes extinct a.s. In the latter case, we provide 
the expectation and the Laplace transform of the absorption time, 
as a functional of the solution of a Riccati differential equation. We 
also show that the quadratic regulatory term allows the LB-process 
to start at infinity, despite the fact that births occur infinitely of- 
ten as the initial state goes to oo. This result can be viewed as an 
extension of the pure-death process starting from infinity associated 
to Kingman's coalescent, when some independent fragmentation is 
added. 

1. Introduction. 

1.1. Presentation of the LB-process. The goal of this paper is to define 
and study the random analogue of a demographic deterministic model which 
is ubiquitous in ecology and widely known as the logistic growth model. It 



Received October 2003; revised September 2004. 

AMS 2000 subject classifications. Primary 60J80; secondary 60J70, 60J85, 92D15, 
92D25, 92D40. 

Key words and phrases. Size-dependent branching process, continuous-state branching 
process, population dynamics, logistic process, density dependence, Ornstein-Uhlenbeck 
type process, Riccati differential equation, fragmentation-coalescence process. 

This is an electronic reprint of the original article published by the 
Institute of Mathematical Statistics in The Annals of Applied Probability, 
2005, Vol. 15, No. 2, 1506-1535. This reprint differs from the original in 
pagination and typographic detail. 



1 



2 



A. LAMBERT 



is defined by the following ordinary differential equation: 

(1) dZ t = bZ t dt - cZt dt, t > 0, 

where b and c are real numbers (c > 0). The quadratic regulatory term has 
a deep ecological meaning as it describes negative interactions between each 
pair of individuals in the population. The stochastic version of the logistic 
process we construct is a density-dependent continuous time branching pro- 
cess, in both continuous and discrete state-space. We call it the "branching 
process with logistic growth" or "logistic branching process," abbreviated 
as LB-process. 

In the discrete state-space setting, individuals give birth as in the pure 
branching case (independently, at constant rate p, to i.i.d. numbers of off- 
spring), but can die either naturally (at constant rate d) or by competition 
pressure [when the total population has size n: at rate c(n — 1), propor- 
tionally to the number of extant conspecifics] . Therefore, when the total 
population has size n, the first birth arrives at rate pn and the first death at 
rate dn + cn(n — 1). This size-dependence complies with the deterministic 
logistic growth model (1), in its formal definition (quadratic death rate), 
as well as in its interpretation (negative interactions between all possible 
pairs). 

In the continuous state-space setting, the LB-process is a Markov process 
with (nonnegative) real values. Its paths are a.s. cdlg, that is, they are 
right-continuous with left-hand limits. The general definition is inspired by 
Lamperti's transform [27] linking continuous-state branching processes and 
Levy processes (i.e., processes with independent and stationary increments), 
and is done by time-changing Ornstein-Uhlenbeck type processes ([34], page 
104). In the special case when the LB-process has a.s. continuous paths, it 
is the unique strong solution (when the initial state is fixed) of a stochastic 
differential equation (SDE) of the following type: 

(2) dZ t = bZ t dt - cZl dt + y^Z~ t dB t , t>0, 

where B is the standard Brownian motion and 7 the so-called Gaussian 
coefficient. In this case, we retain the name coined by Etheridge in [13] where 
she studies it in a spatial setting: "Feller diffusion with logistic growth." 

1.2. Modeling density-dependence. In population biology, the most sim- 
ple process modeling the dynamics of a population is the Malthusian process. 
If Zt £ [0, 00) denotes the total number of individuals at time t, then the 
Malthusian process satisfies dZt = bZ^dt, where b is the mean birth-death 
balance per individual and per time unit. The solutions are straightforward 
exponential functions and when b > 0, they rapidly go to 00, proving useless 
for long-term models. Moreover, this model does not allow populations with 
positive growth to become extinct. 
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This elementary model has a random counterpart, called the branching 
process, where populations may have positive (expected) growth and be- 
come extinct. In their discrete time and discrete state-space form, branch- 
ing processes go back to Lord Francis Galton and Irenee-Jules Bienayme. 
The so-called Bienayme-Galton-Watson (BGW) process is a Markov chain, 
where time steps are the nonoverlapping generations, with individuals be- 
having independently from one another and each giving birth to a (random) 
number of offspring (belonging to the next generation). These (random) off- 
spring all have the same probability distribution. Here, the mean growth is 
geometric, but the process evolves randomly through time, eventually dying 
out or tending to oo, with probability 1. 

Despite its advantage of allowing extinctions while the mean growth is 
positive, the BGW-process shares with the Malthusian process the short- 
coming (from an ecological standpoint) of being able to go to oo. In the de- 
terministic celebrated improvement of the Malthusian process is the 
logistic process (1). It is an elementary combination of geometric growth for 
small population sizes and a quadratic density-dependent regulatory mech- 
anism. 

The main advantage of this model is that Zt converges to a finite limit as 
t — > oo, namely, b/c (if b > 0) or (if b < 0). On the other hand, this model 
does not allow the population to evolve once it has reached its stable state. 
A natural continuation will then be for us to replace geometric growth in 
the logistic equation by random branching (random growth with geometric 
mean). Alternatively, this can be seen as improving the branching process 
by, loosely speaking, adding a quadratic regulatory term to it (and thus 
prevent it from going to oo). 

It is actually a general feature of numerous models in population dynamics 
that the process describing the evolution of the population size through 
time either goes to oo or ultimately dies out [17]. In particular, this is the 
case for any integer- valued Markov chain with as an accessible, absorbing 
state, including (most) BGW processes. As a consequence, the question of 
adding density-dependence to branching processes is not a new concern, and 
answers can roughly be divided into three types (a brief survey in this matter 
can be found in [18]). The first approach sticks to the branching scheme in 
the discrete time setting: at each generation, individuals have i.i.d. random 
numbers of offspring, but their common distribution depends on the current 
population size (see, e.g., [23] and the references therein; see also [14, 35]; a 
special Poisson case is treated in [7] where, despite the title of his paper, the 
author acknowledges that the model is indeed a density-dependent discrete 
time branching process, but has no relation to the logistic growth model (1)). 
The second approach relies on stochastic difference equations, giving rise to 
what is called the controlled BGW process (see, e.g., [16, 21, 33]). The third 
approach consists in generalizing the birth and death rates of the branching 
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process in continuous time in various ways: by considering polynomial rates 
as functions of the population size instead of linear rates, or by setting to zero 
the birth rate outside a compact set of population sizes. This way of modeling 
density dependence is popular among biologists (see, e.g., [11, 29, 30]). 

The present work pertains to the latter approach, but here we are in- 
terested in a specific phenomenon, namely, constant pairwise competition 
pressure, which is a ubiquitous fundamental mechanism in biology in various 
space and time scales, and thus deserves special attention. We will consider 
branching processes in their continuous time form, in both continuous and 
discrete state-space (note, nevertheless, that a similar construction in dis- 
crete time could be done, by sticking to the branching scheme but letting 
the probability of having no child depend on the population size in the same 
fashion) . 

1.3. Statement of results. Because all properties of the integer- valued 
LB-process are common to the real-valued LB-process and not conversely, 
we state the results in the continuous setting. Roughly speaking, these results 
hold under condition (L) that the logarithm of the offspring size at each birth 
time has finite expectation. 

Consider the LB-process with continuous paths presented above (2). We 
will show that this process is a time-change in Lamperti's fashion [27] of the 
Ornstein-Uhlenbeck process, which is itself defined from a positive real num- 
ber c, termed the competititon rate, and the Levy process (^yB t + bt,t>0). 
When the paths need not be continuous, we will consider a general Levy pro- 
cess with no negative jumps instead of the previous Brownian motion with 
drift. The LB-process is thus characterized by c> and a real- valued func- 
tion ip, termed the branching mechanism, which is the Levy exponent of a 
Levy process with no negative jumps. In the case of Feller's diffusion with 
logistic growth (2), ip(X) = ^A 2 - bX, A > 0. 

Let us briefly sum up the results. When the underlying branching mecha- 
nism does not allow the population to decrease in the absence of competition 
pressure, that is, there are no "natural deaths" (ip is then the Levy exponent 
of an increasing Levy process or subordinator), the LB-process is recurrent. 
It is null-recurrent (continuous setting only) if it has zero drift (no "infinites- 
imal births"), Levy measure of finite mass p, further satisfying p < c, and, in 
that case, it converges to in probability. Otherwise it is positive-recurrent 
and its limiting distribution is explicitly given via its Laplace transform 
(Theorem 3.4; Theorem 2.2 in the discrete setting). When there are "nat- 
ural deaths," the LB-process goes extinct a.s. We also display formulae for 
its resolvent measure and its extinction time (via Laplace transform as well, 
Theorem 3.9). 

Let us mention that in the latter case (continuous setting only), the LB- 
process either goes to but remains positive, or is absorbed at in finite 
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time, depending solely on the branching mechanism (Theorem 3.5), that is, 
according to a criterion that does not involve the competition rate c: absorp- 
tion occurs with probability 1 if J°° d\/ip(\) < oo, with probability other- 
wise [where the last inequality has the standard meaning that d\/ip{\) 
exists and is finite for some x]. 

Most of all, we are able to start the LB-process at +oo (Corollary 3.10, 
Theorem 2.3). In the discrete setting this should be compared to the total 
mass process in the standard coalescent starting from infinity (see Section 
2.3). Specifically, we denote by w q the unique positive solution on (0,oo) 
vanishing at +oo, of the following Riccati differential equation: 

II 9 

y - y = -qr , 

where r is some positive function depending on the branching mechanism ip 
and the competition rate c (Lemma 3.8, Lemma 2.1). Namely, set 

0(A) = £dt exp (f*^- ds \ A ^°> 

and (p its inverse function. Then r = (p'/y/ap. The standard LB-process 
(starting from infinity) has entrance law given by 

/ rOM \ 
Eoo(exp(— \Z T )) = exp|^— J w q (s)dsj, \>0,q>0, 

where r is an independent exponential variable with parameter q. In partic- 
ular, the absorption time T a has Laplace transform under P^, 

Eoo(exp(-gT a )) =exp^-^ w q (s)ds\ q>0. 

In addition, its expectation is finite and equal to 

Eoo(T a )= / sr 2 (s)ds. 
Jo 

A shortcoming of our logistic branching process is obviously that it leads 
to much less tractable formulae than in the pure branching case. Though, 
we claim that it is the most natural and realistic model for random self- 
regulatory population dynamics. The simple definition in the discrete set- 
ting can easily be handled by biologists to do simulations or to include these 
dynamics into more complicated models, such as spatial or genetic ones. 
Indeed, in the case when no natural deaths occur, the LB-process survives 
with probability 1 without tending to oo, and might thus provide a fruitful 
framework for randomly evolving large populations, as well as everlasting 
genealogies. Open questions are the following: distributions of independent 
sums of LB-processes, conditioning of the LB-process to survive and quasi- 
stationary distributions, underlying genealogy of the LB-process, distribu- 
tion of coalescence times .... In the pure branching case the first of these 
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questions has an obvious answer, the others have been studied in particular 



The outline of the paper is straightforward. Section 2 is concerned with 
the discrete setting, and is written in a way that is supposed to not scare 
biologists, as it contains the results geared for potential applications. Section 
3 gives a more formal account of the general continuous setting, as well as 
further interesting subtleties that do not appear in the discrete setting (null- 
recurrence, extinction without absorption). Section 4 gathers the proofs of 
various theorems and lemmas. 

2. The discrete setting. In this section we define and study the LB- 
process living on the nonnegative integers. The definition differs from that 
of the LB-process in the continuous setting, but both processes have very 
similar properties. As a consequence, we will state the results in a rather 
straightforward fashion, and proofs will be omitted (we leave the slight adap- 
tations to the reader, and refer her or him to the manuscript, to be found 
on the author's website, for details). 

2.1. A few definitions and a key result. Let us remind the reader that 
in continuous time and in the pure branching case, individuals behave inde- 
pendently as in the BGW-process, but give birth (rate p > 0) and die (rate 
d > 0) at (random exponential) independent times. The number of offspring 
born at each birth time is k with probability irk/ Pi where (iVk,k > 1) is a 
sequence of nonnegative integers such that 



When the population size is n, the total birth rate is thus pn and the natural 
death rate dn. In order to model competition pressure, we then add extra 
deaths at rate cn(n — 1) (c > 0), considering that each of the n individuals 
is in competition with the n — 1 remaining others. These deaths occur as if 
each particle selected another fixed particle at constant rate c and killed it. 
Observe also that these deaths due to competition have the same kernel as 
the pure-death counting process associated to Kingman's coalescent [22], so 
that the LB-process can be viewed as (the mass of ) a combination of coag- 
ulation with independent fragmentation. This viewpoint will be developed 
briefly in Section 2.3. 



in [2, 6, 10, 24, 25, 26, 28, 31]. 




Set 



ip(s) = d-(p + d)s + '^2 TTiS 

i>l 



s G [0, 1] 



which characterizes completely the underlying branching mechanism. 
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The LB-process Z = (Zt,t > 0) associated with tp and the positive real 
number c is thus the (minimal Feller) process with infinitesimal generator 
Q = (qij,i,j > 0), where 



Qij 



iTTj-i, if i > 1 and j > i, 

di + ci(i — 1), if % > 1 and j = i — 1, 

—i(d + p + c(i — l)), if i > 1 and j = i, 

, 0, otherwise. 

From now on, we assume condition (L) is satisfied, that is, 

(L) ^TTilogOO < oo. 

i>l 

Denote by ("Ki)i the tail of the measure (vrj)j, that is, = X^fc 71 "^ & > 1> 
and note that condition (L) is equivalent to J2k(^k/k) < oo. It is easy to see 
that ip(s) = (1 — s)(d — J2k>i^kS k ), which allows, for any s G (0,1] under 
condition (L), the definition 

m ( s > = / ~~T\ T dv - 

Js CV[1-V) 

Standard calculations yield 

(3) exp(m( S )) = p S ~ d / c exp ( £ ^ Y s £ (0, 1], 



where 



P = exp - 5^ 



fe>i 



Define the nonnegative decreasing function 0: (0,1] i— ► [0,£) by 

0( s ) = [ e m ^dv, w€ (0,1], 

where £ = Jq 1 e m ^'-' € (0, +oo]. As a consequence of (3), we can assert that 

£ = oo d> c. 

The mapping is a bijection, whose inverse on [0, £) will be denoted by 
In particular, 

(p'(s) = -exp(-moip(s)), sG[0,(). 
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Lemma 2.1. Assume d^O, and denote by (e' q ) the following Riccati 
differential equation: 



«,) y' -y 2 = -qr 2 

where 



r(.) = -_!$2l s£(0>0 . 

Vc^(s)(l -y>(a)) 

For any positive q, (e' q ) has a unique nonnegative solution w q defined on 
(0,£) and vanishing at £— . In addition, w q is positive on (0,£), and for 
any s sufficiently small or large, w q (s) < yfqr{s). As a consequence, Jq w q 
converges, and w q decreases initially and ultimately. 

When the coordinate process starts from x, its law will be denoted by P x . 

2.2. Results. The behavior of the LB-process depends on whether d is 
positive or zero. As previously, assume condition (L) holds. 

When d = 0, recall (3) to see that the function exp(m) can be extended 
continuously to [0, 1] . It is not difficult to show that one can define the 
probability measure v on N by 

exp(m(s)) = ^ v i sl ~ X i sG [0,1]. 

i>l 

Theorem 2.2. When d = 0, the LB-process is positive-recurrent in N, 
and converges in distribution to the probability measure ji defined by 

\j>l 



In the binary-splitting case (p = TT\), this limiting distribution is that of a 
Poisson variable of parameter p/c conditioned on being positive 

Pi \- e -p/ c i\ ' 

Recall the function w q defined in Lemma 2.1. Let T a denote the absorption 
(extinction) time. 

Theorem 2.3. When d>0, the LB-process goes extinct a.s. Moreover, 
the probabilities (¥ x ,x > 0) converge weakly, as x — > oo, to the law of 
the standard LB-process, or logistic branching process starting from infinity. 
Under P^, the extinction time T a has Laplace transform 

Eoo(exp(-gr a )) =exp^-^ w q {z)dz^j, q>0, 
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its expectation is finite and equal to 




2.3. The link with fragmentation-coalescence processes. Recall that when 
ip = and c> 0, the LB-process is merely the pure-death process associated 
to Kingman's coalescent [22], where the coagulation rate per pair of ob- 
jects is a constant equal to 2c. As a consequence, the LB-process can be 
regarded in general as the mass of a coalescent process with independent 
fragmentation (corresponding to the branching mechanism): components of 
the partition coalesce at rate 2c, split independently into a random number 
k > 2 of subcomponents at rate itk—ii an d may also disappear spontaneously 
at rate d (erosion). 

It is known that the pure-death process (p = 0) can be started at infinity, 
which stems from the fact that arrival times of coalescence events (coagula- 
tion from i objects into i — 1 objects) have summable expectations, propor- 
tional to — 1)) • Then Theorem 2.3 asserts that it is still true even after 
adding fragmentation. Nevertheless, under P^, fragmentation events occur 
infinitely often with probability 1 as t \ 0. Indeed, the following probabili- 
ties are not summable: 

Pi (a birth occurs before a death) = — - — ; ~ — as % — > oo. 

n ; (p + d)i + ci(i- 1) ci 

Let us mention a few recent results in this vein. 

In the discrete reaction-diffusion model studied in [3], particles perform 
independent random walks, die spontaneously and locally coalesce or split 
with the same kernels as ours (in the special case of binary branching). 

In [4], all exchangeable coalescence-fragmentation (EFC) processes living 
in the partitions of N are characterized. In a case of quick coalescence and 
slow fragmentation (Kingman's coefficient is nonzero and fragments always 
come in finite numbers), the mass of the EFC process has the same transition 
kernel as our (discrete) logistic branching process. This allows Berestycki 
to use Theorem 2.3 to show that these EFC processes "come down from 
infinity," in the sense that the number of components of the partition is 
finite at any positive time t a.s., even though the process starts from dust (at 
t = 0, it is equal to the infinite collection of singletons). Note, however, that 
Theorem 2.3 was stated under the assumption that d ^ 0. To be convinced 
that the EFC process still comes down from infinity even when d = 0, observe 
that when there are n > 2 objects, the total death rate is cn{n — 1) > |n(n — 
1) + |n, which is the total death rate associated to the LB-process with 
competition rate c/2 and death rate c/2. 

Other works have considered size-biased fragmentation-coalescence pro- 
cesses (on the partitions of an interval; rates are proportional to the lengths 
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of the subintervals), when fragmentation is binary. These partitions have 
infinitely many components and a ubiquitous equilibrium measure in this 
setup is the Poisson-Dirichlet distribution [9, 32]. 

2.4. Convergence to the logistic Feller diffusion. A convergence result 
is easily obtained that links the discrete and continuous settings. Namely, 
Feller's diffusion with logistic growth is the limit of a sequence (indexed 
by n £ N) of discrete logistic branching processes. Rigorously, denote by 

(Z[ n \t >0) the Markov process living on n X N, started at n 1 \nx\ , stopped 

(n) (n) 

at 0, and whose transition kernels are as follows. The integer iV f = nZ\ 
is a binary-splitting {tti = p) LB-process with parameters (indexed by n) 
Pn = ^n 2 + Xn, d n = ^n 2 + 5n, c n = c, where c, 5, 7, A, are positive constants. 
In other words: 

1. z[ nS> is incremented by 1/n at rate + X)n 2 z[ n \ 

2. z[ n ^ is decremented by 1/n at rate (^n + 5)n 2 z\ n \ 

3. Z^ is decremented by 1/n at rate cn 2 Z^ n \z^ — n" 1 ). 

Let 6 = A — 5. Then standard results [20] show that as n — > 00, the sequence 
{z[ n \t > 0) n converges weakly to the Feller diffusion with logistic growth 
started at x and solution of (2), 

dZ t = bZ t dt - cZl dt + ^/^Z~ t dB t , t > 0. 
3. The continuous setting. 

3.1. Preliminaries. In the 1950s [19], the analogue of the BGW-process 
was defined in continuous time and continuous state-space (CB-process). 
All CB-processes are cdlg, and their jumps (if any) are a.s. positive. Those 
whose paths are continuous satisfy an SDE of the following type: 

dZ t = bZ t dt + y/jZ~ t dB t , t > 0, 

so that the definition of the LB-process in this case is clearly (2). When 6 = 0, 
the dynamics are referred to as quadratic branching. When 6 = and 7 = 4, 
Z is the celebrated Feller diffusion or squared Bessel process of dimension 
0. 

Next we remind the reader of a celebrated result of Lamperti [27] which 
relates CB-processes and Levy processes with no negative jumps. This result 
will allow us to give the appropriate definition of the branching process with 
logistic growth in a second section. 

A Levy process is a cdlg Markov process with independent and stationary 
increments. A Levy process with continuous paths is the sum of a (scaled) 



THE LOGISTIC BRANCHING PROCESS 



11 



Brownian motion and a deterministic drift (possibly zero). Let X be a real- 
valued Levy process with no negative jumps. Let To be the first hitting time 
of zero by X. Then define 

,*AT ds 

Vt= t>o, 

Jo -&s 

and (Ct,t > 0) its right-inverse. Lamperti's result then states that if 

Z t = X(C t ), t>0, 
then Z is a CB-process, and 

rt 



C t = I Z s ds, t>0. 
Jo 



Conversely, any CB-process Z is a time-changed Levy process: if C is defined 
as above, and rj is the right-inverse of C, then Z o 77 is a Levy process with 
no negative jumps killed when it hits 0. 

In the special case of quadratic branching, Lamperti's result predicts that 
if we time-change a Feller diffusion (i.e., a quadratic CB-process) in this 
fashion, we obtain (up to a multiplicative constant) a killed Brownian mo- 
tion. It is indeed easy to check that if Z% = x + Jq \fZ~ s dB s , then Z o rj is a 
local martingale with increasing process Jq* Z s ds = t, and thus is a standard 
Brownian motion. Using this last argument, it is elementary to prove the 
following proposition concerning the LB-process. 

Proposition 3.1. As previously, assume Z is a Feller diffusion with 
logistic growth, that is, a diffusion solving (2), 

dZ t = bZ t dt - cZf dt + y/^Z~ t dB t , t > 0. 

Define rj as the right-inverse of C, where C% = f^Zgds. Then the process 
R = Z on solves the following SDE: 

dR t = dX t - cR t dt, t > 0, 

where X is a Brownian motion with drift, namely, ^^^(Xt — bt,t> 0) is 
the standard Brownian motion. 

Conversely, let R be the (strong) solution of the last SDE where X is a 
Brownian motion with a possible drift, Tq its first hitting time of and C 
the right-inverse of n, where 

f tAT ds 

Vt= t>0. 
Jo R s 

Then Z = X o C is a diffusion process killed when it hits 0, and solves an 
SDE of type (2). 



This last proposition will enable us to extend the definition of branching 
process with logistic growth to any kind of branching mechanism. 



12 



A. LAMBERT 



3.2. Construction of the logistic branching process. Let X denote a spec- 
trally positive Levy process (i.e., with no negative jumps). The branching 
mechanism function that interests us is the Levy exponent ip of X defined 
by 

E(exp(-AX t )) = exp(^(A)), t,A>0. 
It is specified by the Levy-Khinchin formula ([5], Chapter VII) 

i>{\) =a\ + J A 2 + /°°(e" Ar - 1 + Xrl r<1 )U(dr), A > 0, 
2 Jo 

where a is a real number, 7 a positive real number termed the Gaussian 
coefficient of X, and LT is a positive cr-finite measure on (0, 00) such that 
Jo°(l A r 2 )H(dr) < 00, termed the Levy measure of X. 
We now provide the definition of the LB-process. 

definition 3.2. For any positive real number c and any spectrally pos- 
itive Levy process X with Levy exponent tjj, we define the logistic branching 
process LB(^,c) associated with tp and c starting from x > as follows. 

Let R denote the unique strong solution, starting from x, of the following 
SDE: 

(4) alR t = dX t - cR t dt, t > 0, 

To its first hitting time of and C the right-inverse of rj, where 

f tAT ds 

Vt= TTi *>0- 

JO Ms 

Then the LB-process Z is the Feller process defined by 

R(C t ), if < t < t/oo, 
0, if 7700 < 00 and t > r]^. 



Z t 



If A (resp. Q) denotes the infinitesimal generator of X (resp. Z), then, for 
any differentiable / in the domain of A, 

Qf(z) = zAf(z)-cz 2 f'(z), z>0. 

Before continuing further, we prove some of the facts stated in this last 
definition. SDEs of type (4) are well studied, see, for example, [34], pages 
104-113. They have strong solutions which are cdlg homogeneous strong 
Markov processes, known as Ornstein-Uhlenbeck type processes. An explicit 
formula for the unique strong solution of (4) with initial state x is 

(5) R t = xe~ ct + exp(-ci) / exp(cs) dX s , t > 0. 

Jo 
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By standard theory of Markov processes (see, e.g., [12]), Z is then a cdlg 
time homogeneous strong Markov process. Furthermore, it is easily seen that 
the infinitesimal generator U of R is given, for any differentiable / in the 
domain of A, by 

Uf(z)=Af(z)-czf'(z), zeR. 
As a consequence, for any time t and initial condition x > 0, with s = C u , 

K x Qf Qf{Z u ) du\ = E x Uf(Z u )Z u du^ 

= K x (£ t Uf(R s )ds^ 

= K x (f(R(C t ))) - f(x)=E x (f(Z t )) - f{x), 
which shows, indeed, that Q is the infinitesimal generator of Z . 

3.3. Properties of the logistic branching process. In this section we are 
interested in the law of the LB-process, and particularly in its long-term 
behavior. This behavior depends on whether AC is a subordinator or not. We 
recall that a Levy process with positive jumps is called a subordinator if it 
has increasing paths a.s. In that case, its paths have finite variation, its drift 
coefficient is nonnegative, and its Gaussian coefficient is zero. 

In contrast to the pure branching case, LB-processes eventually go to 
with probability or 1. In reference to population biology, this probability is 
called the extinction probability (otherwise, absorption probability). When 
X is a subordinator, this probability is 0, and the LB-process is recurrent. 
When X is not a subordinator, it is 1, and we characterize the law of the ex- 
tinction time. Unless otherwise specified, proofs of statements of this section 
are postponed to Section 4. 

From now on, we assume condition (L) is satisfied, that is, 

(L) E(log(X 1 ))<oo. 

We recall that condition (L) is equivalent to f°° \og(r)H(dr) < oo [5]. Ac- 
cordingly, it is easily seen that one can define 

m(X)= [ X ^ds, A>0. 
Jo cs 

In the next statements, we consider the case when X is a subordinator. We 
then denote by S > its drift coefficient, so that ([5], Chapter III) 

POO 

ij)(X) = -SX- II(dr)(l- e~ Xr ), A > 0. 
J o 

We also introduce its Levy tail S, that is, n(y) = Il(dr), y > 0. 
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Lemma 3.3. Assume X is a subordinator satisfying (L). Then m can 
be expressed as 

-m(A) = ~A+/ (l-e~ Xr )^dr, A > 0, 
c Jo cr 

and the following equation 

u(dr)e- Xr = exp(m(A)), A>0 



defines a unique probability measure v on (0, oo). This probability measure 
is infinitely divisible and v((5/c, oo)) = 1. 

Remark. Note that v is also self-decomposable. An r.v. S is said to be 
self-decomposable if for any a < 1, there is an r.v. S a independent of S such 
that 

aS + S a == S. 

In particular, any self-decomposable distribution is infinitely divisible. Con- 
versely, a real infinitely divisible distribution is self-decomposable if its Levy 
measure is of type \x\~ 1 k(x) dx, where k is increasing on (— oo,0) and de- 
creasing on (0, oo). 

We introduce condition (d), where p is defined as 



poo 

/ n(dr) < oo. 
Jo 



We thus say that (d) holds iff (at least) one of the following holds: 

• 6^0, 

• p = oo, 

• c < p < oo. 

Theorem 3.4. Assume X is a subordinator satisfying (L). 

(i) Assume (d). Then the probability v defined in Lemma 3.3 has / °° r _1 x 
v(dr) < oo. The process LB(tp,c) is positive-recurrent in (5/c, oo) and con- 
verges in distribution to the probability measure \i whose size-biased distri- 
bution is v, that is, 



fi(dr) = ( / s 1 u(ds) ) r l v(dr), 

\J(S/c,oo) J 



r > 0. 



' (<5/e, oo) 

In particular, the expectation of the stationary probability is 



L 



rfi(dr) = ( / s l v(ds) ) < oo. 

(5/c,oo) \J(5/c,oo) 
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(ii) Assume (<9) does not hold. Then the process LB(tp,c) is null-recurrent 
in (0, oo) and converges to in probability. 

From now on, X is assumed not to be a subordinator. In the next theorem, 
we claim that the extinction probability is 1. We set a criterion to establish 
whether the process remains positive or is absorbed, that is, reaches in 
finite time. Note that the criterion for absorption does not depend on c and 
is the same as for the branching process (c = 0) [15]. 

Theorem 3.5. Assume X satisfies (L) and is not a subordinator. Then 
the process LB(ip,c) goes to a.s., and if T a denotes the absorption time 

T a = mf{t>0:Z t = 0}, 

then P(T a < oo) = 1 or according to whether f°° dA/ip(X) converges or 
diverges. 

The last two theorems will be proved in Section 4. 

Next we intend to give deeper insight into the law of the LB-process. We 
will simultaneously be able to define the standard LB-process or LB-process 
starting from infinity. 

Define 9 : [0, +oo) i-> [0, +oo) by 

0(A) = [ X e m ^ds, A>0. 
Jo 

Recall [5] that when X is not a subordinator, liminf^oo X~ 1 ip(X) > 0. As a 
consequence, the mapping 9 is a strictly increasing bijection, whose inverse 
will be denoted by 99. In particular, 

ip (A) = exp(— m o (p(X)), A > 0. 

Define also for any nonnegative x and positive the Laplace transform Gq^x 
of the g-resolvent measure of the LB-process (Z t ,t > 0) started at x, 

POO 

G q , x {\)= dte- qt E x (e- XZt ), A > 0. 
Jo 

Note that from Definition 3.2, if A still denotes the infinitesimal generator of 
the Levy process X, then for any differentiable real function in the domain 
of A, 

roo 

q / dte- qt E x (f(Zt)) 
(6) J ° 

roc 

= f{x) + / dte-^ x {Z t Af[Z t )-cZlf{Z t )). 
Jo 
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Before stating the main theorem of this section, we display three lemmas, the 
last of which is the key result for the theorem. The first two help understand 
where the third one comes from. These lemmas are also proved in Section 
4 (except Lemma 3.8 whose proof is replaced by the more general proof of 
Lemma 2.1). 

Lemma 3.6. As a function of X £ (0,oo), G q>x is twice continuously 
differentiable and solves the second-order linear differential equation (E q ), 
where y is the unknown function and X is the scalar variable 

(E q ) -cXy" + ^y' + qy = e- x \ 

Lemma 3.7. 1} 
ated to (Eg), then: 



Lemma 3.7. If (Eq ) is the homogeneous differential equation associ- 



(i) For any solution (I, f q ) of (E q ), for any open subinterval Iq of I on 
which f q does not vanish, the function g q defined as 

solves on Iq the Riccati differential equation (e q ) 

^ - y+y+ 7x y = -\- 

(ii) For any solution (J,g q ) of (e q ), the function h q defined as 

h q (\) = e- mo ^g q o ip (\), AG#(J) 
solves on 6 (J) the Riccati differential equation (e' q ) 

( I \ I 2 2 

(4) y - y = -or , 

where 

LEMMA 3.8. For any positive q, the Riccati differential equation (e' ) 
has a unique nonnegative solution w q defined on (0, oo) and vanishing at 
oo. In addition, w q is positive on (0, oo), and for any X sufficiently small or 
large, w q (X) < y/qr(X). As a consequence, J w q converges, and w q decreases 
initially and ultimately. 



Now we are able to state the main results of this section. 
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Theorem 3.9. Recall the function w q defined in Lemma 3.8. For any 
nonnegative A and positive q, an expression for G q ^ x {\) is given by 



MX) 

qG q>x (X) = l- dte 

(7) 



0(A) r e(A) 
x I" ds qr 2 {s){\ - e- x ^)e~ ft w « {z) dz . 



In particular, if J°° d\/ip(\) converges, then the absorption time T a is a.s. 
finite and has Laplace transform under ¥ x , 

E x (exp(-qT a )) 

i_ r dte -fr w ^ dz 



(8) 

dsqr 2 (s)(l - e - x ^)e-It w " iz)dz , q > 0. 
In addition, its expectation is finite and is equal to 

roc 

E x (T a )= / dssr 2 (s)(l-e- x ^) 

(9) J ° 

= / ^ (1 _ e -**) e -m(t) / dse m(s)_ 
Jo ct Jo 

Corollary 3.10. The probabilities (¥ x ,x > 0) converge weakly, as x — > oo, 
to the law Poo of the standard LB-process or logistic branching process start- 
ing from infinity. The standard LB-process is the Markov process that has 
the same transition kernels as the LB-process considered previously, and en- 
trance law given by 

/ ,0(A) \ 

(10) E 00 (exp(-AZ T )) = exp( -J w q {s)ds), \>0,q>0, 

where r is an independent exponential variable with parameter q. In partic- 
ular, if f°° d\/ip(\) converges, then under P M the absorption time T a is a.s. 
finite and has Laplace transform 

(11) E 00 (exp(-qT a )) = exp^- J w q (s)ds^j, q>0. 
In addition, its expectation is finite and equal to 

roo roo A+ rt 

(12) E OQ (T a ) = sr 2 (s)ds = — e~ m ^ / dse m( - s \ 

Jo Jo ct Jo 
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Proof of Theorem 3.9. As in Lemma 3.7, define g q as 

g q (s) = e m ^w q o9(s), s>0. 

Then the function g q solves the Riccati differential equation (e 9 ) and by 
Lemma 3.8 is integrable at 0, for 

A r0(\) 

g q (s)ds= / w q (t)dt, A > 0, 



by the change t = 9(s). Then define the functions F q x and K q ^ x as 
F q>x (X) = G q>x (X)exp(J^ g q (s)ds^, A > 0, 

K q7X (X) = F^ x (X)exp- (m(X) + 2j\ q (s)d s y A > 0. 
By Lemma 3.6 and after some algebra, we get 

cXK q x {\) = — exp — ^Ax + m(A) + J g q (s)ds^j, A > 0. 
Next integrate the last equation to obtain 



K q>x (X)=K + [ 

J A 



M e -tx-m(t)-f*g q (s)ds 
A ct 

oo ,J2i 



dv ^) e -x V {v)~f w q {u)du^ A>Q) 

Je{x) c<p(v) 

where the last equality stems from the changes u = 9(s) and v = 0(t), and 
Kq is an integration constant. Therefore, by the same kind of change, we 
get 

Jo(X) c<p(s) 

Check that this last function is integrable at 0. Indeed, F qx (X) ~ K qiX (X) 
as A \ 0, and the first expression we got for K qtX provides the equivalent 
K q ^ x {\) ~ c ln(A). Then because F qjX (0) = G qjX (0) = 1/q, the integration 
of F' reads 



1 r°W o r* 

F q>x (X) = -+K dte 2 Jo 
q Jo 



= 2 Jo w i^ dz 



+ / dteJo w ^ z ' az I dgTL-fLe-xrw-JtW", A > 0. 

ap(s) 
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We thus get the following expression for G q ^ x : 

, . , r 0(\) , „ , fd(X) „ r t , . , 

G/ \ \ — 1 — / w„(z)dz , T y —I w„(z)dz I j, 2 I w„(z)dz 
qx{A) = q e J o qy ' + Koe J o qy ' / ate J o qy ' 

JO 

•J w q (z)dz dte J o w q (z)dz dsr 2^ e -x V (s)-J t w q (z)dz^ 



+ e 

o 

A> 0. 

Recall that qr 2 = w q — w' q , so that an integration by parts yields 

f(?(A) ft , , , /-oo r s 

dte J Q w q (z)dz j dsqr 2^ e -J t Wq (z)dz 

= / dtw q (t)e J o qy ' 
Jo 

= e J w q {z)dz_ l _ 

Since gG^o is constant equal to 1, the constant integration Kq must be 0. 
In addition, the previous display yields for any nonnegative A, 

rem 



qG q , x (\) = e-Io Wq ^ dz 



r0(X) f0(\) 
+ / dte-Jt Wq{z)dz 
o 



(14) x I" dsqr 2 (s)e- x ^e-fi Wq{z)dz 



0(A) f 9(A) 



x / dsqr 2 (s)(l-e- x ^)e-J t Wq{z)dz , 



which ends the proof of (7). For (8), recall that T a is a.s. finite upon the 
condition given in Theorem 3.5. Observe that 



fOO 

(A)=/ dte~^ x {l Zt=0 + e~ XZt l Zt>Q ), A>0. 
Jo 



By a double application of Lebesgue's convergence theorem, this quantity 
converges to 

p oo 

/ dte- <lt F x (Z t = 0), 
Jo 

so that, integrating by parts the last quantity, 

lim qGq JX) = E x (exp(-qT a )). 

A— >oo 
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For the expectation of T a (9), we have to handle the Laplace exponent H (t, A) 
of the marginal distribution of the LB-process. Fubini's theorem and Beppo 
Levi's theorem, indeed, yield 

;>oo ;>oo 

lim / (l-H(t,X))dt = E l Zt ^ dt = E{T a ). 

A^ooJo JO 

It is clear (see proof of Lemma 3.6) that 

dH ,,^dH ,d 2 H 

^ = -^(A)^ + cA^, t,X>0, 

which yields, by integrating ^e -m ( A ) w.r.t. A, 

^ = -e m(A) r-^(M)e-^ ; t,X>0. 
oX J\ cs at 

Integrate this last equation w.r.t. A and then w.r.t. t, and recall that H(t, 0) = 
1, to get, by Fubini's theorem, 

/*oo rX roo 

/ (l-H(t,X))dt= / due m(u) / — e - m( - s \l-e- sx ), A > 0. 

J JO Ju CS 

The result follows by letting A go to oo and applying Lebesgue's convergence 
theorem. An elementary variable change yields (9). □ 

Proof of Corollary 3.10. Equations (10) and (11) of the corollary 
follow from the same arguments that supported the theorem. 

Let us prove the convergence of E x (exp(— XZ T )) to (12). Observe that if r 
is an independent exponential r.v. with parameter q, then ~E x (exp(—XZ T )) = 
qGq )X {X). The convergence then stems from (14) and the integration by parts 
used in (13). Namely, 

< E,.(exp(-AZ T )) - e ~C X)w ^ dz 



< 



0(A) rB(^) , \ , /*oo r s , . , 

dte-Jt e -Mt) j dsqr 2 (s)e-Jt w * {z)dz 

0(A) ,em 



dte-ft w ^ dz e- x ^w q (t). 
Now the function 



is integrable on [0, 0(A)], so that, by Lebesgue's convergence theorem, the 
last displayed upper bound converges to as x — > oo . 

Finally, the existence of (and weak convergence to) the LB-process start- 
ing from infinity stem from the previous convergence of Laplace transforms 
and standard theory [8]. □ 
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4. Proofs. 

4.1. Proofs of Lemma 3.3, Theorems 3.4 and 3.5. 
Proof of Lemma 3.3. First observe that 

fA /7 C / roc 



-m(X)= [ —(5s+rn{dr)(l-e- rs ) 
Jo cs V Jo 



5X /*oo / i \ /" Ar , 1 - e" 



o 



h / n(dr) / du 

c Jo Jo cu 

= Vf(l-e--)2M*. 
c Jo cr 

Note that the positive measure A defined by A(dr) = (cr)~ 1 Tl(r) dr, r > 0, 
is, under condition (L), a Levy measure. Indeed, 

roc 

A(dr) = c" 1 / ln(u)U(du) < oo, 
Jl 

poo 

rk(dr) = c" 1 / (1 A u)U{du) < oo, 



o Jo 

so that one can define Y as the subordinator with Laplace exponent — m ([5], 
Chapter III). As a consequence, the probability measure v is well defined 
as the law of Y\ (so it is infinitely divisible), and the measure A always has 
infinite mass so that Y\ > 5/c a.s. □ 

Proof of Theorem 3.4. To see that 

POO 

(15) (d)& / r~ l u{dr) < oo, 

Jo 

we use the following equality (where both sides can be infinite): 

) POO 

r - 1 v(dr)= / e m( - x U\, 



so that (15) reduces to (d) f°° dX exp(m(A)) < oo. 

Recall that m(A) = f (cs)~ 1 ^(s) ds < 0, A > 0. It is clear from Lemma 
3.3 that whenever 5 ^ 0, f°° dX exp(m(A)) is finite. 

Next assume that IT has finite mass p. Since the case when II has infinite 
mass can easily be derived from the case p > c by a truncation argument, 
it only remains to prove that when 5 = 0, f°° dX exp(m(A)) < oo iff p> c. 
Then recall Lemma 3.3 and write, for any A > 1, 



roo 

-m(A)= / (1 - e" Ar )^ dr. 
Jo cr 
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Next pick x > 0, divide (0,oo) into (0,x/X\, (x/X,x] and (x,oo), and change 
variables (u = Ar), to get 

— m(A) — / — — dr 
Jx/x cr 

cu 

w ' du + (1 - e~ Ar )— ^-dr. 



P n(r) 

lim — m(X) — / dr 



'x cu J x cr 

Since II is positive and decreases from ft(0+) = p, an appeal to Lebesgue's 
convergence theorem yields 

" x n(r) 

Ix/X 

(16) =p r (1 _ e -» ) ^_,r e -»^ + r5M tir 

Jo cu J x cu cr 

= k(x), 

and /c(x) is finite. Now the same properties of IT ensure that 

(17) HMw A )< T ^ldr< ^-hi(X), x>0,X>l. 

c " Jx/x cr c 

Using (16) and the left-hand side inequality in (17), we get 
limsupe m(A) A n ^/ c < exp(-fc(x)) < oo. 

A— >oo 

Assume first that p > c and pick x small enough to have ti{x)/c > 1, then 
the last equation shows that f°° dX exp(m(A)) < oo. 

Using (16) and the right-hand side inequality in (17), we get 

liminf e m(A) A p/c > exp(-k{x)) > 0. 

A— >oo 

If P < c, the last equation then shows that J°° dX exp(m(A)) = oo. This ends 
the proof of (15). 

We next point out that, (d) holding or not, Z is recurrent in (5/c, oo). 
Indeed, it is known [34] that R oscillates in (<5/c, oo) and converges in dis- 
tribution to z/, so that To = oo a.s., lim^oo = oo a.s., and C is defined on 
the entire [0, oo), so Z = RoC also oscillates in (S/c, oo). 

Let us characterize the invariant measures. According to Definition 3.2, 
the infinitesimal generator Q of Z satisfies, for any positive A, 

Qe\(z) = zAe\(z) - cz 2 e' x {z) = {ip(X) + cXz)ze\(z), z>0, 
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where e\(z) = e Xz . Now for any positive measure ( on [0, oo), let x denote 
the Laplace transform of rQ(dr). Then Q is invariant for Z iff £Q = and 

(Q = (Qe x = VAX) 

C(dz)(ip(X) + cXz)ze x (z) = VA>0 

^(A)x(A)-cA X / (A) = 0, A>0. 

The resolution of this ordinary differential equation proves that x(A) = 
/cexp(m(A)) for some constant k, so that 

(18) r((dr) = ku(dr), r > 0, 

with a possible Dirac mass at for £. We can now prove the dichotomy 
between (i) (convergence in distribution) and (ii) (convergence to in prob- 
ability) . 

(i) Thanks to (15), since (d) holds, we know that r~ l v(dr) < oo. 
Equation (18) then ensures that \i is the unique invariant probability mea- 
sure. 

(ii) In the case when (<9) does not hold, Jq°° r~ 1 u{dr) = oo and all invariant 
positive measures on (0, oo) are nonintegrable at 0+, so that limt^oo F(Zt > e) 
for any positive e. 

Recall that Z oscillates in (0, oo) to conclude it is null-recurrent. □ 



Proof of Theorem 3.5. In the case when X is not a subordina- 
tor, the Ornstein-Uhlenbeck type process R still oscillates [34] with sta- 
tionary measure u, but this time v charges the negative half-line so that 
lim inft—xx) Rt < a.s., and To < oo a.s. It is thus clear that lim^oo Z% = 
a.s., the absorption depending on whether the integral Jq° ds/R s (= T a ) con- 
verges or not. To check the criterion for absorption given in the theorem, we 
make use of the function G q)X . First define the positive function g qx as 

%x(A) = -^(A), A>0. 

^q,x 

In the rest of the proof, we will drop the indices whenever it is not ambiguous 
to do so. Recall that 

E(exp(-gr a )) = lim qG(\). 

A— >oo 

Since qG(X) = exp(— Jq g(s) ds), it suffices to show that 

/*oo /*oo 

/ g(X)dX and / dX/i^(X) have the same nature. 
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Thanks to Lemma 3.6, it is not difficult to prove that g solves 

cX(-y / + y 2 ) + iPy = q-—-, A > 0. 

We first show that G decreases more slowly than any negative exponential 
function of parameter, say, z 

e~ zX 
lim — — ■ = 0. 

A^oo G(X) 

To this aim, define 

a = sup{i > : Z t > z/2}, 
which is finite a.s. since Z goes to 0. Next observe that 

G(X) > / dt e~ qt E(e" AZt ,t>a) 
Jo 

> e ~Xz/2 



poo 

/ dte- qt F{t >a) 
Jo 



= e~ Xz/2 q- l ¥.{(T qa ), 

which proves our claim. Since X is not a subordinator, X/i/j(X) is bounded 
from above (and positive for sufficiently large A), so that 

-Ax 

ft(A) =^<wr 0(e 1 asA ^°°' 

We will use this last comparison along with the following rearrangement of 
the differential equation that g solves 

(19) i>(\)(g + h)(X)-q-cXg'(X) = -cXg 2 (X), A > 0. 

On the vector space Cp of continuous real functions defined on [0, oo) with 
(at most) polynomial growth at oo, define the scalar product (•,•) by 

(u,v) = ^ ^ \ v ; -, u,veC P . 

In particular, if 1 denotes the constant function equal to 1 and I the identity 
function on [0, oo), then (1,1) = 1 and 

by the Cauchy-Schwarz inequality. Thus, the positive function g is decreas- 
ing on [0, oo). If I > is its limit at oo, since limsup^^ g'(X) = and 
limA^oo X/tp(X) £ [0, oo), (19) yields / < 0, so that / = 0. 
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Thanks to this same equation, and because I = 0, note that 
lim g(xy\g(X)+h(X)-^(Xr 1 (q + cXg'(X))) = - lim -£-g(X)=0. 

A^oo A-+00 ^(A) 

As a consequence, 

/oo 
g(X) dX and 

/oo 
(V'(A) _1 ((/ + cA</(A)) — h(X)) dX have the same nature. 

But h(X) = o(e~ Ax '/ 2 ), liniA^oo \/i/>(\) G [0,oo) and lim A ^oc I 5(A) = Z = 0, 
so that 

/•OO /"OO 

/ A^(A) _1 </(A) dA and / h(X)dX both converge, 
which reduces the comparison to J°° dX/ip(X). □ 
4.2. Proofs of Lemmas 3.6, 3.7 and 2.1. 

Proof of Lemma 3.6. Recall (6), and apply it to the function e\ 
defined as e\(z) = e~ Xz . It is well known that e\ is in the domain of A and 
that, for any positive A, 

Ae x (z)=i;(X)ex(z), z>0. 

As a consequence, 

;>oo 

qG q , x {X) = e x (x) + / dte- qt EMX)Z t e- XZt +cXZ?e~ XZt ) 
Jo 

= e\(x) - ip(X)G' q x (X) + cXGq X (X), 
which concludes the proof. □ 

PROOF of Lemma 3.7. Since f q solves the homogeneous equation (E q h ^) 
associated to (E q ), we have 

-cXf' q '{X) + V>(A)^(A) + qf q (X) = 0, AG/. 

It is an old trick to consider g q = -f' q /f q , for then f' q '/f = -g' q + (f' q /f q ) 2 - 
For A G Iq [f q {X) 7^ 0], divide the previously displayed equality by —cXf q (X), 
and get 

-g' q (X) + g 2 q (X) + ^g q W = ± X G J„. 
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For any g q satisfying the last equation on some open interval J, we consider 
the change of function h g (X) = e- mo ^g q o <p(\), A G 9 (J), so that 

h' q (X) = e-^W (_ m ' o <p(\)g q o </>(A) + 5 ; o 99(A)) 

= e - 2m o ¥ ,(A) ^2 QV9(A) 



K{\) 



cip(\) 

which indeed proves that h q solves (e') on 8 (J). □ 



Proof of Lemma 2.1. 



Notation. For any real number t and real functions f,g which do not 
vanish in a pointed neighborhood of t, we will write f(s) ~ g(s) [resp./(s) x 
g(s)] as s — ► i to mean 

lim = 1 (resp. = some nonzero, finite real number). 

s->t g( S ) 

Recall that (e^) is the following Riccati differential equation: 

(e' g ) y' -y 2 = -qr 2 , 

where 

and £ < 00 <£4> d < c. By the Cauchy-Lipschitz theorem (see, e.g., Corollaries 
1 and 2, Chapter 7.2, page 93 in [1]), for any T £ (0,£) and a£R, there is 
a unique maximal solution (I, /) to (e^) with boundary condition f(T) = a. 
We denote by f(r,a) this maximal solution and by I(r,a) the corresponding 
maximal interval. We will use repeatedly, and sometimes implicitly, the fact 
that the orbits of different solutions never intersect, therefore inducing a 
total ordering on any set of solutions defined on the same interval. 



Let us give the outline of the proof. We state three lemmas and then 
provide the proof of Lemma 2.1. The first lemma studies the graph of r and 
shows, in particular, that r is monotone on neighborhoods of 0+ and £— , 
and that lim i _ > o+ r (t) — +00. The second one is a technical result ensuring, 
in particular, that, for any T 6 (0, £), the maximal interval I(t,o) is the entire 
(0,£). The third lemma asserts that the family (/(t,o))t is uniformly increas- 
ing and bounded from above, allowing us to define w q as their increasing 
limit as T — > £ — . 
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The rest of the section will then be devoted to proving that w q solves 
(e' q ), and that it is the unique nonnegative solution vanishing at £— , and to 
studying its behavior at 0+ and £— . 

Lemma 4.1. The function r is monotone on a neighborhood of £— and 
monotone increasing on a neighborhood o/0+. Furthermore, 

lim r(s) = +oo, 



lim r(s) = 




ifd< c/2, 
ifd = c/2, 
ifd>c/2. 



We point out that when d > c, one has £ = +oo, lim s _> +00 r(s) = 0, and 
lim s _>o+ r (s) = +oo, which is the same situation as in the continuous setting. 
The sequel of the proof in the case when d > c is thus exactly the same in 
the continuous setting. 

Proof of Lemma 4.1. Set 

-2m(s) 

H(s)=r 2 oe(s) = — se(0,l), 

cs(l — s) 

and recall (3) to see that 

#0) ~ -w>s 2d/c ~ 1 ass^0+ and H (s) ~ . 1 . ass->l-. 
cp^ c(l — s) 

Then differentiate -ff to check that it is monotone on neighborhoods of 
and £, so that the same holds for r, by monotonicity of 9. □ 

LEMMA 4.2. (i) For any T £ (0,£) and a£l, 

a < =>• [T, £) C ^(t,o) anc ^ /(T.a) ^ s negative on [T, £), 

a > =>■ (0,T] C ^(T,a) an d f(T,a) ^ s positive on (0, T]. 

(ii) If d < c/2, there is a neighborhood V of £ s-uc/t £/iai /or any Tg V 
and aGl, 

a<T/qr(T) ==>■ [T, £) C -fyr.a) anc ^ nm /(T,a)(^) exists and is finite. 

Proof. First note that for any T £ (0,£), f\ T ,o)( T ) = and f[ Tfi) (T) = 
—qr 2 (T) < 0. Consequently, a maximal solution of (e' g ) vanishes, at most, 
once on its interval of definition. 
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(i) Let a < 0. We just saw why / = f(T,a) remains negative for t > T. 
Now assume that the supremum a of I(r,a) IS m Since / is a maximal 
solution and is negative on (T, £), we deduce that lim^o-. f(t) = — oo. But 
this would mean that / is below —^fqr on some neighborhood of a, and thus 
increasing. This shows that a = £. The same arguments hold for the other 
implication (a > 0). 

(ii) Assume a G [— y/qr(T), y/qr(T)]. When d < c/2 (Lemma 4.1), r in- 
creases on some neighborhood V of £ < oo, and goes to oo at £— . For this 
reason, / = /(T, a ) decreases (/' = f 2 — qr 2 ), so it remains below ^/qr, but 
never reaches —^/qr, since at a hitting time of —^fqr, f would vanish, 
whereas —^fqr decreases. So / is defined, as claimed, at least on [T,£) even 
if a > 0, and for any t in a neighborhood of £, / 2 (t) < qr 2 (t). The same holds 
when a < — y/qr(T), for then / increases until it reaches —yfqr. 

To get the existence of a finite limit at £— , it is thus equivalent, by /' — 
/ 2 = — qr 2 and |/'| < qr 2 , to prove that f^r 2 < oo. To see this, use (3) and 
get, as s — > £— , 

~ (f _ s )V(l-d/c) SQ that r 2 (s) x (e _ s) -(l-2i/c)/(l-d/c)_ 

The result follows by observing that if d G (0, c/2), then ^ ( — ^'^)- 

□ 

As a consequence of the previous lemma, for any T £ (0, £), /(t,o) is defined 
on the entire (0, £), that is, I(t,o) = (0)0- Next set 

«T = /(T,0)> TG(0,O- 

Lemma 4.3. Tne sequence (kt)t is uniformly increasing and uniformly 
bounded from above on (0, £) . One can i/ien define the function w q as the 
uniformly increasing limit 

w a = lim t kt- 

The function w q is positive on (0,£) and lim^g_ w q (t) = 0. 

Proof. The orbits (kt)t are totally ordered because they never inter- 
sect, and increase as T increases because k' t (T) = —qr 2 (T) < 0. 

We next intend to prove that the family (kt)t is uniformly bounded from 
above. 

(i) d> c. In this case, £ = oo, and there is a to such that r decreases 
on [to,oo). Thanks to Lemma 4.2(i), ho = f(t lhy /qr(t )) 1S defined (at least) 
on (0,to] an d is positive. Since different orbits cannot intersect, for any 
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Tg (0,oo), kt is below ho on (0,to]- Conclude by noticing that on [to, 00), 
K T < \/Q r (otherwise kt would increase on [to, 00), which contradicts Lemma 
4.2(0). 

(ii) d£ [c/2,c). Recall Lemma 4.1 to see that r can be extended contin- 
uously to (0,£]. In the same vein as in Lemma 4.2(i), the maximal solution 
^ of (e') vanishing at £ is defined on (0, £]. It is then obvious that kt is 
below Kg for any T < £. 

(iii) d < c/2. Thanks to Lemma 4.2(h), there is a to close enough from 
£ such that /io = f(t ,y/qr(to)) has a positive limit at £— . This function /i is 
defined on (0,£) and since it never vanishes, it never intersects any kt- As 
a consequence, the family (kt)t is bounded from above by /to- 

Now for any t G (0,£) and T > t, Kr{t) > [Lemma 4.2 (i)] . Since u> g (t) 
is the increasing limit of kt(£) as T increases, w q {t) > 0. The function w q 
vanishes at £— by continuity of the flow. □ 

We hereafter conclude the proof of Lemma 2.1. We first prove that for 
any positive q, w q is the unique nonnegative solution of (e'„ ) defined on (0,£) 
and vanishing at £— . 

Let us show the uniqueness of such a solution. Assume h and k are two 
nonnegative solutions of (e' ) on a neighborhood of £, who both vanish at £— . 
If there is some to such that h(to) 7^ A; (to), say /i(to) > fc(to), then u = h — k 
is positive on (to,£) and so is w = h + A;. Since h and A; are both solutions of 
(e' q ), u and v are such that u' = uv, so that u is positive and increasing on 
(to,£). This contradicts the fact that lirnt_*g_ u(t) = 0. 

To show that w q is the solution of (e' q ), fix e in (0, £), and set 77 = w g (e) and 
h e = /( £ ,^)- Next suppose that for some to £ -f( £ ,r/), *o > e and h e (to) / w q (to). 
Since is by construction the supremum of a family of orbits indexed 
by a continuous set, it is necessary that /i £ (to) > w q (to). Pick then some 
a G (w q (to),h £ (to)). We know from Lemma 4.2(i) that / = /(f 0)O ) is defined 
at least on (0,to]. Since /(to) > w q (to), /(to) > ^r(^o) f° r all T G (0, £). 
This last inequality certainly holds for all t G (0, to], and thus, in particular, 
f(e) > w q (e) = rj. Reasoning in the same fashion with h £ , we see that /(e) < 
h £ (e). But h £ (e) = 77, which brings the contradiction (77 < /(e) < 77), and 
proves that u> g (t) = h £ (t) for all t G ^( e ,rj), t > e. Now by the same reasoning 
as in the proof of Lemma 4.2, h £ is defined (at least) on [e, £), since otherwise 
it would go to +00 in finite time and there would, indeed, be some to such 
that h e (to) > w q (to). As a consequence, w q = h £ on [e,£). This ends the 
proof, since e is arbitrarily small. 

We complete the proof of Lemma 2.1 by studying the behavior of w q at 
0+ and 

First note that because w' q = w 1 q — qr 2 , w q increases when it is above y^r, 
and that w' q vanishes when both graphs intersect. Recall r is positive and 
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monotone on neighborhoods of and £. As a consequence, the graph of w q 
may intersect that of ■ s fqr at most once on each of these two neighborhoods. 
On the other hand, it cannot remain above yjqr on a whole neighborhood: 

(i) of 0+, because it would be increasing but tend to oo as t \ 0+; 

(ii) of £— , because it would be increasing and tend to as t f £— . 

Now since w q < ^fqr on neighborhoods of 0+ and £— , it decreases on these 
neighborhoods, and /jio, converges because Jq r does so. Indeed, 

r(s) ~ —FT ^ =, ^ ass^0+ and r(s) ~ — |L=^= ass^(-, 
Vc(l - y(s)) V C( A S ) 

which are both integrable on their respective neighborhoods, since ip has 
finite limits at 0+ and £— . □ 
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